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Abstract 

The structural, electronic and magnetic properties of small Fe^Rh^ clusters having A'' = m+n < 
8 atoms are studied in the framework of a generalized-gradient approximation to density-functional 
theory. The correlation between structure, chemical order, and magnetic behavior is analyzed as a 
function of size and composition. For N = m + n <Q a thorough sampling of all cluster topologies 
has been performed, while for N = 7 and 8 only a few representative topologies are considered. 
In all cases the entire concentration range is systematically investigated. All the clusters show 
ferromagnetic-like order in the optimized structures. As a result, the average magnetic moment 
per atom JIj^ increases monotonously, which is almost linear over a wide range of concentration 
with Fe content. A remarkable enhancement of the local Fe moments beyond 3 //s is observed 
as result of Rh doping. This is a consequence of the increase in the number of Fe d holes, due 
to charge transfer from Fe to Rh, combined with the extremely reduced local coordination. The 
Rh local moments, which are important already in the pure clusters {N < 8) are not significantly 
enhanced by Fe doping. However, the overall stability of magnetism, as measured by the energy 
gained upon spin polarization, increases when Rh is replaced by Fe. The composition dependence 
of the electronic structure and the influence of spin-orbit interactions on the cluster stability are 
discussed. 

PACS numbers: 75.75.-|-a, 36.40.Cg, 75.50.Bb, 73.22.-f 



I. INTRODUCTION 



Alloying elements with complementary qualities in order to tailor their physical behavior 
for specific technological purposes has been a major route in material development since the 
antiquity. Cluster research is no exception to this trend. After decades of systematic studies 
of the size and structural dependence of the most wide variety of properties of monoelement 
particles, the interest has actually been moving progressively over the past years towards 
investigations on finite-size binary alloys.'^ The magnetism of transition-metal (TM) clusters 
opens numerous possibilities and challenges in this context. '2^3 Yot example, one would 
like to understand how to modify the magnetic characteristics of clusters, in particular the 
saturation magnetization and the magnetic anisotropy energy (MAE), as it has been done 
in solids. This would indeed allow one to design new nanostructured materials from a 
microscopic perspective. Nevertheless, it also true that controlling composition, system size, 
and magnetic behavior sets serious difficulties for both experiment and theory. 

Pure TM clusters such as Fe^, Con and Ni^r show spin moments, orbital moments, and 
MAEs that are enhanced with respect to the corresponding periodic solids.'^^'^ Still, the 
possibilities of optimizing the cluster magnetic behavior by simply tuning the system size 
have been rather disappointing, particularly concerning the MAE, which remains relatively 
small — despite being orders of magnitude larger than in solidd^^' — due to the rather weak 
spin-orbit (SO) coupling in the 3d atoms. This is one of the motivations for alloying 3d 
TMs with 4d and 5c? elements which, being heavier, are subject to stronger SO interactions. 
In this context it is useful to recall that large nanoparticles and three dimensional solids of 
these elements are non-magnetic. However, at very small sizes the 4d and 5d clusters often 
develop a finite spontaneous low-temperature magnetization, due to the reduction of local 
coordination and the resulting d-band narrowing.f^^^ The first experimental observation of 
this important finite-size effect has been made by Cox et al. by performing Stern- Gerlach- 
deflection measurements on RhTv clusters. In this work the average magnetic moments per 
atom Jljq = 0.15-0.80yUB have been experimentally determined for < 30-50 atoms.'^^ In 
view of these contrasting features one expects that 3d-4:d and 3d-5d alloy clusters should 
show very interesting structural, electronic and magnetic behaviors. 

The purpose of this paper is to investigate the ground-state properties of the small FeRh 
clusters in the framework of Hohenberg-Kohn-Sham's density functional theory.l^Sl Besides 
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the general interest of the problem from the perspective of Sd-Ad nanomagnetism, these 
clusters are particularly appealing because of the remarkable phase diagram of FeRh bulk 
alloys.!^ In the case of FesoRhso the magnetic order at normal pressure and low temperatures 
is antiferromagnetic (AF). As the temperature increases this a" phase undergoes a first order 
transition to a ferromagnetic (FM) state, the a' phase, which is accompanied by a change 
in lattice parameter. The corresponding transition temperature T"'"" increases rapidly 
with increasing external pressure P, eventually displacing the FM a' phase completely for 
P > 7 GPa (T"'"" ~ 290K for FesoRhso at normal pressure). Moreover, T"'"" decreases 
very rapidly with decreasing Rh content. At low pressures the FM a' phase undergoes a FM 
to paramagnetic (PM) transition at (Tc — 670K).^^ In addition, the properties of a-FeRh 
bulk alloys have been the subject of first principles and model theoretical investigations.'^ 
In particular these show that the relative stability of the FM and AF solutions depends 
strongly on the interatomic distances. Such remarkable condensed-matter effects enhance 
the appeal of small FeRh particles as specific example of 3d-4d nanoscale alloy. Investigations 
of their magnetic properties as a function of size, composition, and structure are therefore 
of fundamental importance. 

The remainder of the paper is organized as follows. In Sec. [IT] the main details of the 
theoretical background and computational procedure are presented. This includes in par- 
ticular a description of the strategy used for exploring the cluster energy landscape as a 
function of geometrical conformation and chemical order. The results of our calculations for 



FeRh clusters having < 8 atoms are reported in Sees. Ill and IV, First, we focus on the 
interplay between structure, chemical order and magnetism in the most stable geometries 
for different cluster sizes. Second, we analyze the concentration dependence of the cohesive 
energy, the local and average magnetic moments, and the spin-polarized electronic structure. 
Finally, we conclude in Sec. |V]with a summary of the main trends and an outlook to future 
extensions. 



II. COMPUTATIONAL ASPECTS 

The calculations reported in this work have been performed in the framework of Hohenberg- 
Kohn-Sham's density functional theory,^*^ as implemented in the Vienna ab initio simulation 
package ( VASP) .I^S The exchange and correlation energy is described by using both the spin- 
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polarized local density approximation (LDA) and Perdew and Wang's generalized-gradient 
approximation (GGA).'231 The VASP solves the spin-polarized Kohn-Sham equations in an 
augmented plane-wave basis set, taking into account the core electrons within the projector 
augmented wave (PAW) methodP This is an efficient frozen-core all-electron approach 
which allows to incorporate the proper nodes of the Kohn-Sham orbitals in the core region 
and the resulting effects on the electronic structure, total energy and interatomic forces. 
The 4s and 3d orbitals of Fe, and the 5s and 4d orbitals of Rh are treated as valence 
states. The wave functions are expanded in a plane wave basis set with the kinetic en- 
ergy cut-off Emax = 268 cV. In order to improve the convergence of the solution of the 
selfconsistent KS equations the discrete energy levels are broadened by using a Gaussian 
smearing a = 0.02 eV. The validity of the present choice of computational parameters has 
been verified.'^Sl The PAW sphere radii for Fe and Rh are 1.302 A and 1.402 A, respectively. 
A simple cubic supercell is considered with the usual periodic boundary conditions. The 
linear size of the cell is a = 10-22 A, so that any pair of images of the clusters are well 
separated and the interaction between them is negligible. Since we are interested in finite 
systems, the reciprocal space summations are restricted to the P point. 

Although the potential advantages of alloying magnetic 3d elements with highly-polarizable 
4d or 5d elements can be grasped straightforwardly, the problem involves a number of serious 
practical challenges. Different growth or synthesis conditions can lead to different chemical 
orders, which can be governed not just by energetic reasons but by kinetic processes as well. 
For instance, one may have to deal with segregated clusters having either a Ad core and a 
3d outer shell or vice versa. Post-synthesis manipulations can induce different degrees of 
intermixing, including for example surface diffusion or disordered alloys. Moreover, the inter 
atomic distances are also expected to depend strongly on size and composition. Typical 
TM-cluster bond-lengths are in fact 10-20% smaller than in the corresponding bulk crystals. 
Taking into account that itinerant 3d-electron magnetism is most sensitive to the local and 
chemical environments of the atoms,'^^'^^^ it is clear that controlling the distribution of the 
elements within the cluster is crucial for understanding magnetic nanoalloys. 

Systematic theoretical studies of binary-metal clusters are hindered by the diversity of 
geometrical conformations, ordered and disorder arrangements, as well as segregation ten- 
dencies that have to be taken into account. This poses a serious challenge to both ffist- 
principles and model approaches. In order to determine the interplay between cluster struc- 
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ture, chemical order and magnetism in FeRh clusters we have performed a comprehensive 
set of electronic calculations for clusters having N < 8 atoms. In the present paper we focus 
on the most stable cluster structure and magnetic configuration, which are determined by 
exploring the ground-state energy landscape.'^ This is a formidable task, since one needs 
to consider a large, most possibly complete and unbiased set of initial structures. Such a 
thorough geometry optimization must include not only the representative cluster geometries 
or topologies, but also all relevant chemical orders. This requires taking into account all 
distributions of the Fe and Rh atoms for any given size and composition. These two aspects 
of the problem of determining the structure of nanoalloys are discussed in more detail in the 
following. 

The different cluster topologies are sampled by generating all possible graphs for N < 6 
atoms as described in Ref. (SSj (see also Ref. I^Tjl . For each graph or adjacency matrix it is 
important to verify that it can be represented by a true structure in D < 3 dimensions. 
A graph is acceptable as a cluster structure, only if a set of atomic coordinates Ri with 
i = 1, . . . , N exists, such that the interatomic distances Rij satisfy the conditions Rij = Ro 
if the sites i and j are connected in the graph (i.e., if the adjacency matrix element Aij = 1) 
and Rij > Rq otherwise (i.e., if Aij = 0). Here Rq refers to the nearest neighbor (NN) 
distance, which at this stage can be regarded as the unit of length, assuming for simplicity 
that it is the same for all clusters. Notice that for < 4 all graphs are possible cluster 
structures. For example, for N = 4, the different structures are the tetrahedron, rhombus, 
square, star, triangular racket and linear chain.'^Sl However, for > 5 there are graphs, 
i.e., topologies, which cannot be realized in practice. For instance, it is not possible to 
have five atoms being NNs from each other in a three dimensional space. Consequently, for 
A^ > 5 there are less real structures than mathematical graphs. The total number of graphs 
(structures) is 21 (20), 112 (104), and 853 (647) for A^ = 5, 6, and 7, respectively.lSSl 

For clusters having A^ < 6 atoms all these topologies have indeed been taken as starting 
points of our structural relaxations. Out of this large number of different initial configura- 
tions the unconstrained relaxations using VASP lead to only a few geometries, which can 
be regarded as stable or metastable isomers. For larger clusters (A^ = 7 and 8) we do not 
aim at performing a full global optimization. Our purpose here is to explore the interplay 
between magnetism and chemical order as a function of composition for a few topologies that 
are representative of open and close-packed structures. Taking into account our results for 
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smaller sizes, and the available information on the structure of pure FeN and RhAr clusters, 
we have restricted the set of starting topologies for the unconstrained relaxation of FeRh 
heptamers and octamers to the following: bicapped trigonal bipyramid, capped octahedra, 
and pentagonal bipyramid for N = 7, and tricapped trigonal bipyramid, bicapped octahe- 
dra, capped pentagonal bipyramid and cube for = 8. Although, the choice of topologies 
for = 7 and 8 is quite restricted, it includes compact as well as more open structures. 
Therefore, it is expected to shed light on the dependence of the magnetic properties on the 
chemical order and composition. 

The dependence on concentration is investigated systematically for each topology of 
FcmRhn by varying m and for each size A^ = m + n < 8, including the pure Fe^ and 
Rh^v limits. Moreover, we take into account all possible non-equivalent distributions of the 
m Fe and n Rh atoms within the cluster. In this way, any a priori assumption on the 
chemical order is avoided. Obviously, such an exhaustive combinatorial search increasingly 
complicates the computational task as we increase the cluster size, and as we move away 
from pure clusters towards alloys with equal concentrations. Finally, in order to perform the 
actual density-functional calculations we set for simplicity all NN distances in the starting 
cluster geometry equal to the Fe bulk valu^ Rq = 2.48 A. Subsequently, a fully uncon- 
strained geometry optimization is performed from first principles by using the VASP.E2lThe 
atomic positions are fully relaxed by means of conjugate gradient or quasi-Newtonian meth- 
ods, without imposing any symmetry constraints, until all the force components are smaller 
than the threshold 5 meV/A. The convergence criteria are set to 10^^ eV/A for the energy 
gradient, and 5 x 10~^ A for the atomic displacements.!^ The same procedure applies to 
all considered clusters regardless of composition, chemical order, or total magnetic moment. 
Notice that the diversity of geometrical structures and atomic arrangements often yields 
many local minima on the ground-state energy surface, which complicates significantly the 
location of the lowest-energy configuration. 

Lattice structure and magnetic behavior are intimately related in TMs, particularly in 
weak ferromagnets such as Fe and its alloys.^^ On the one side, the optimum structure and 
chemical order depend on the actual magnetic state of the cluster as given by the average 
magnetic moment per atom JIj^ and the magnetic order. On the other side, the magnetic 
behavior is known to be different for different structures and concentrations. Therefore, in 
order to rigorously determine the ground-state magnetic properties of FeRh clusters, we have 
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varied systematically the value of the total spin polarization of the cluster 5*2 by performing 
fixed spin-moment (FSM) calculations in the whole physically relevant range. Let us recall 
that Sz = (z^t — '^4.)/2 where ^'^(^'4.) represents the number of electrons in the majority 
(minority) states. In practice we start from the non-magnetic state (S"™" = 0) and increase 
5*2 until the local spin moments are fully saturated, i.e., until the Fe moments in the PAW 
sphere reach /ifg — and the Rh moments //Rh — 2.5 fiB (typically, S^'^^ > 3N/2). The 
above described global geometry optimizations are performed independently for each value 
of Sz- These FSM study provides a wealth of information on the isomerization energies, the 
spin-excitation energies, and their interplay. These are particularly interesting for a subtle 
magnetic alloy such as FeRh, and would therefore deserve to be analyzed in some more detail. 
In the present paper we shall focus on the ground-state properties by determining for each 
considered FcmRh^ the most stable structural and magnetic configuration corresponding to 
energy minimum as a function of and of the atomic positions.^Sl 

Once the optimization with respect to structural and magnetic degrees of freedom is 
achieved, we derive the binding energy per atom Eb = [mElFe) +nii^(Rh) — E{Fem.^h.n)]/N 
in the usual way by referring the total energy E to the corresponding energy of m Fe and 
n Rh isolated atoms. Moreover, for each stationary point of the total energy surface (i.e., 
for each relaxed structure having a nearly vanishing IIV-EH) we determine the vibrational 
frequencies from the diagonalization of the dynamical matrix. The latter is calculated from 
finite differences of the analytic gradients of the total energy. In this way we can rule 
out saddle points to which the local optimization procedure happens to converge on some 
occasions. Only configurations which correspond to true minima are discussed in the follow- 
ing. Finally, a number of electronic and magnetic properties — for example, the magnetic 
energy AE^ = E{Sz = 0) — E{Sz), the local magnetic moments /Xj integrated within the 
Wigner-Seitz (WS) or Bader atomic cells of atom i^SSEsl g^^^j ^j-^g gpj^ polarized density of 
electronic states (DOS) Pcr(e) — are derived from the self-consistent spin-polarized density 
and Kohn-Sham spectrum. 

III. STRUCTURE AND MAGNETISM 

In this section we discuss the ground-state structure, chemical order, binding energy, and 
magnetic moments of Fe^Rhn clusters having N = m + n < 8 atoms. The main empha- 
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sis is here on understanding how the various electronic, structural and magnetic properties 
depend on the chemical composition of the alloy. First, each cluster size is analyzed sep- 
arately, since a strong dependence on is expected in the small size, non-scalable regime. 
Comparisons between the various are stressed by means of cross-references between dif- 
ferent subsections. In addition the main trends as a function of size and concentration are 
summarized in Sec. IIVI 

A. FeRh dimers 

Despite being the simplest possible systems, dimers allow to infer very useful trends on 
the relative strength, charge transfers and magnetic order in the various types of bonds 
which are found in FeRh alloy clusters. The results summarized in Table [T] show that the 
FeRh bond yields the highest cohesive energy, followed by the Rh2 bond, the Fe2 bond being 
the weakest. The particular strength of the heterogeneous bond is confirmed by the fact that 
the corresponding vibrational frequency is the highest. The bond length, however, follows 
the trend of the atomic radius which, being larger for Rh, gives (iRhRh > c^FeRh > dpePe- 
Quantitatively, the binding energy per atom E^*^^ = 1.35 eV obtained for Fe2 within the 
GGA is smaller than the LDA result E^^^ = 2.25 eVSSl although it still remains larger than 
the experimental value E^"^^ = 0.65 eV reported in Ref. EU The calculated vibrational 
frequency i/o(Fe2) = 288 cm~^ is consistent with previous experimental results [i^o(Fe2) = 
299.6 cm~^ from Ref. El] and z/o(Fe2) = 300 ± 15 cm^^ from Ref. |52]. Our result for Eb 
and JJjy of Rh2 coincide with previous GGA calculations by B. V. Reddy et al^ These are 
however larger than the experimental values -E^^'^*(Rh2) = 1.46 eV derived from Knudsen 
effusion,'5SlE5'P*(Rh2) = 0.70 ±0.15 eV derived from resonance Raman in Ar matrice^ and 
£'^''''*(Rh2) = 1.203 eV derived from the resonant two-photon ionization.'^ The calculated 
vibrational frequency of z/o(Rh2)*~^^^ = 224 cm~^ should be compared with the experimental 
value uo(Rh2Y''P^ = 283.9 cm"^ reported in Ref. EH 

The stability of magnetism, as measured by the difference in the total energy AEm of 
the non-magnetic {Sz = 0) and optimal magnetic solutions, is largest for Fe2 and smallest 
for Rh2. The same trend holds for the average magnetic moment per atom which decreases 
linearly from JI2 = to 2^b as one goes from Fe2, to FeRh, to Rh2. These average 
magnetic moments per atom correspond to a full polarization of all d electrons in the WS 
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TABLE I. (Color online) Structural, electronic and magnetic properties of FeRh dimers. Results 
are given for the binding energy Eb (in eV), the magnetic stabilization energy ^E^ji — E(^Sz — 
0) — E{Sz) (in eV), the average interatomic distance da/B (in A) between atoms a and /3 (a, /3 = Fe 
or Rh), the average spin moment per atom JZjy = "^Sz/N (in ^b)-, the local spin moment (in 
/is) at the Fe or Rh atoms, and the vibrational frequency uq (in cm~^). 
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spheres: v^t — ^ for Fe and i/^ ~ 8 for Rh, where Ud stands for the number of valence d 
electrons of the corresponding atom. The local magnetic moments = Fe and Rh), 

are obtained by integrating the spin density within the PAW spheres which have the radius 
'"PAw(Fe) = 1.3 A for Fe and rpAw(R-h) = 1.4 A for Rh. In the pure dimers, the local 
moments /ipe = 2.82/iB and /^Rh = 1.83/iB are close to the respective total moment per atom 
7i2 = and 2^b-, which indicates that the spin-density m(r) = n-^{f) — n^if) is quite 
localized around the atoms. Actually, the differences between and /Ijv give a measure of 
the small spill-off effect in m{r). Taking this into account, the results for fia in the FeRh 
dimer seem quite remarkable. Here the Fe local moment is significantly enhanced with 
respect to the Fe2 or Fe-atom value, while the Rh moment is reduced by a similar amount 
(A/xpe = 0.52 fiB and A/jRh = -0.50 hb, see Table [l]). This is mainly the consequence of 
a transfer of d electrons from Fe to Rh, which allows the Fe atom to develop a larger spin 
moment, due to the larger number of d holes. This occurs at the expense of the moment at 
the Rh atom, which has less d holes to polarize. An integration of the electronic density in 
the Bader cell^ shows that 0.33 electrons are transferred from the Fe to the Rh atom in 
FeRh. This behavior is qualitatively in agreement with the higher Pauling electronegativity 
X of the Rh atom (xFe = 1-83 and XRh = 2.28)^ 
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B. FeRh trimers 



The results for trimers are summarized in Table |TTj As expected, the lowest energy isomers 
are found to be triangles for all compositions. According to our calculations the ground 
state of Rha is an equilateral triangle (Dsh) with Eb = 2.31 eV, bond length d = 2.37 A 
and average magnetic moment JI^ = IfiB- The local magnetic moments fia = 0.93/2^ in 
the WS cells align parallel to each other and are almost as large as JI^. These results are 
consistent with those reported in previous GGA studies of Rhs {Eb = 2.35 eV, d = 2.45 A 
and 7I3 = l/is).'^ A single Fe substitution yields an isosceles FeRh2 with an elongated base 
composed of the two Rh atoms. Notice that the bond-length c/pthRh = 2.57 A is larger than in 
Rhs. The linear isomer of the form Rh-Fe-Rh, i.e., with only FeRh bonds, lies 0.4 eV above 
the optimal structure. It is the only true local minimum among the linear FeRh trimers. The 
other linear structures (Rh-Rh-Fe, Fe-Rh-Fe, and Fe-Fe-Rh) are all found to be saddle points 
connecting triangular minima of the potential energy surface (PES). Further Fe substitution 
yields a isosceles Fe2Rh in which the FeFe bond is the shortest. One observes, as in the 
dimers, that the interatomic distances follow the trends in the atomic radii. Finally, for 
Fe3, the calculated lowest-energy structure is a Jahn- Teller distorted isosceles triangle with 
two longer bonds {di2 = di^ = 2.30 A) and a shorter one (^23 = 2.07 A). The calculated 
average magnetic moment of Fes is T^s = 3.33/xb- These results coincide with previous GGA 
studied^ predicting du = dis = 2.33 A and 6/23 = 2.09 A. In contrast, LDA calculation^ 
yield an equilateral Fes with average magnetic moment Jl^ = 2. 66/15 and d = 2.10 A. 
By using the spin-polarized LDA, we also obtain an equilateral triangle similar to the one 
reported in Ref. In contrast, in the GGA one finds that the equilateral triangle is 
unstable with respect to a Jahn- Teller distortion. The isosceles shape of Fes can therefore be 
interpreted as a consequence of exchange and correlation effects. Moreover, we have analyzed 
the GGA Kohn-Sham spectrum in the equilateral structure and found a high degeneracy at 
the Fermi energy, which is consistent with the interpretation that the distortion is triggered 
by a Jahn- Teller effect. 

Concerning the composition dependence of Eb one observes a non-monotonous behavior 
as for N = 2, which indicates that the FeRh bonds are the strongest. The lowest vibrational 
frequency follows a similar trend, despite the larger mass of Rh. Notice that FeRh2 is 
somewhat more stable than Fe2Rh, since the bonds between Rh atoms are in general stronger 
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TABLE II. (Color online) Structural, electronic and magnetic properties of FeRh trimers. Results 
are given for the binding energy per atom Eb (in eV), the magnetic stabilization energy per atom 
AEm = [E{Sz = 0) — £'(5z)]/N (in eV), the average interatomic distance dajS (in A) ordered from 
top to bottom as dpeFe, c^FeRh and dRhRh, the average spin moment per atom = 2Sz/N (in fis) 
, the local spin moment /Uq (in hb) at the Fe or Rh atoms, and the lowest vibrational frequency 
uq (in cm^^). 
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than between Fe atoms. Finally, one may also notice that the energy gain AE^ associated 
to magnetism only plays a quantitative role in the relative stability of triangular and linear 
FeRh2. AEm is actually larger for the linear chain than for the triangle. Therefore, the later 
remains the most stable structure even in the non-magnetic case, although with somewhat 
different bond lengths. 

The average magnetic moment per atom JI^ amounts to IfiB for Rhs- In the alloys it 
increases monotonously with Fe doping, reaching Jt^ = 10/3/iB for Fes- The local mag- 
netic moments fia always show a FM-like coupling. They are all identical in Rhs, which 
is consistent with the Cs point-group symmetry. In the pure clusters /i^ is always close to 
7I3. This indicates that the spin polarization is dominated by electrons occupying localized 
states and that spill-off contributions are not important. For example, in the case of Fes, 
one finds /ii = 3.23/iB and /i2 = yUs = 2.87yUB, the latter corresponding to the pair of atoms 
forming the shorter bond. On the other side, the average local moments /iRh = 0.93/iB 
in Rhs should be compared with 7Z(Rhs) = IfJ^B- As soon as FeRh bonds are present, for 
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mixed compositions, the local Fe moments are enhanced beyond 3/1^- This is mainly due 
to a charge transfer from Fe to Rh, leading to an increase in the number of Fe d holes as 
already observed in the dimer. Quantitatively, the local /ipe and /iRh in mixed trimers are 
similar, though somewhat smaller, to the corresponding values in the FeRh dimer. Notice, 
moreover, the enhancement of the Rh local moments in Fe2Rh and FeRh2 as compared to 
pure Rhs. This reflects the importance of the proximity of Fe on the magnetic behavior of 
the Rh atoms. 



FeRh tetramers 



The most stable FeRh tetramers are all tetrahedra and the first low-lying isomers are 



rhombi (see Table III). The distribution of the atoms within the optimal topology does not 
play a role since all sites are equivalent in a tetrahedron. In the case of Rh4 we obtain a 
nonmagnetic undistorted tetrahedron having Eb = 2.75 eV and bond length d = 2.45 A. 
The closest isomer is found to be a bent rhombus with an average bond length d = 2.35 A. 
Similar results have been obtained in previous studies on Rh clusters.l^ Notice, however, 
that Bae et alW^ have obtained a bend rhombus as the ground-state structure for Rh4 also 
by using VASP. This discrepancy is likely to be a consequence of the different choice of the 
pseudopotential and cutoff energy E^ax- In our calculations we considered the PAW method 
and Emax = 268 eV, while in Ref. l62lone used ultrasoft pseudopotentials and Emax = 205.5 
eV. 

The binding energy of the alloys shows a characteristic non-monotonous dependence on 
concentration, which was also found in smaller clusters. In fact Fe2Rh2 and FeRha are the 
most stable tetramers with Eb = 2.74 eV and Eb = 2.76 eV, respectively. This confirms 
that the FeRh bonds are the strongest. It is worth noting that these trends are not altered 
qualitatively if magnetism is neglected, i.e., if one considers Eb for Sz = 0. In addition, 
it is interesting to follow how Eb changes from Rh4 to Fe4. The stability of the clusters 
can be qualitatively related to the number of homogeneous and heterogeneous bonds by 



counting them for each of the clusters shown in Table III For instance, FeRhs, which is 
the most stable composition, has 3 FeRh and 3 RhRh bonds. Replacing a Rh by an Fe 
to obtain Fe2Rh2 implies replacing 2 RhRh bonds by a stronger FeRh and a weaker FeFe 
bond. Therefore, Eb does not change significantly. The fact that Eb depends weekly on 
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TABLE III. (Color online) Structural, electronic and magnetic properties of FeRh tetramers as in 
Table IE 
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composition for Rh rich tetramers shows that FeRh and RhRh bonds are comparably strong 
in these clusters. 

Concerning the magnetic moments one observes a approximately linear dependence of 
as a function of Fe content. In general, the substitution of a Rh by and Fe atom results 
in an increase of the total moment 25*2 by 3 or 4/iB, or equivalently, A/Z^ = (0.75-l)/iB (see 



Table III). The magnetic order is always FM-like. In the alloys the local moments /ipe show 
the above mentioned enhancement, which is due to a Fe-to-Rh c?-electron charge transfer 
that increases the number of d holes and allows for the development of /ipe — 3.2-3.4/1^. In 
addition, the presence of Fe in Fe^Rhn enhances the Rh local moments as compared to pure 
Rh4. 



D. FeRh pentamers 



In Table [IV] the results for FeRh pentamers are summarized. Although all possible cluster 
topologies (20 structures) were considered as starting geometries for each composition, only 
the most highly coordinated trigonal bipyramid (TBP) and the square pyramid (SP) are 
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found to be most stable geometries. The low coordinated structures transform into compact 
structures after the relaxation. Except for Rhs, which optimal structure is a SP, all the other 
FeRh pentamers have the TBP as ground-state geometry. The trend in the composition de- 
pendence of the binding energy Eb of pentamers confirms the behavior we started to observe 
for = 4. Indeed, in the Fe-rich limit Eb increases rapidly with increasing Rh content, as 
the weakest FeFe bonds are replaced by FeRh bonds. Later on, near 50% concentration and 
in the Rh-rich limit, the composition dependence is weak since FeRh and RhRh bonds are 



comparably strong (see Table IV). In particular for Rh-rich compositions, replacing Fe by 
Rh atoms no longer results in weaker binding. In other words, FeRh bonds are no longer 
primarily preferred. This is possibly a consequence of the increasing coordination number, 
which enhances the role of electron delocalization and band formation, thus favoring the 
larger Rh hybridizations. 

The calculated optimal structure of Rhs, a square pyramid, coincides qualitatively with 
previous DFT calculations.'^ Nevertheless, we obtain a binding energy that is 0.07 eV per 
atom lower than in Ref. (53] Substituting one Rh atom by Fe yields FeRh4 and changes the 
optimal cluster topology to the more compact TBP. The SP remains a local minimum of 
the ground-state energy surface, which is only 3 meV per atom less stable than the optimal 
TBP geometry. The average magnetic moment 7Z(FeRh4) = 1.2^b is enhanced with respect 
to Rhs due to the contribution of a large Fe local moment /ipe = 3.31/iB. Notice that the 
Rh moments are no longer enhanced as in the smaller FeRhAr_i but significantly reduced: 
A^Rh = 0.62/iB for the apex atoms and //Rh = 0.52/2^ for the Rh atoms sharing a triangle with 
the Fe. This is of course related to the fact that the ground-state Sz is relatively low. The 
effect is even stronger in the case of the SP isomer of FeRh4. Here we find two Rh moments 
/^Rh = 0.43/iB that couple parallel to the Fe moment, one very small Rh local moment 
AiRh = 0.05/is, and an antiparallel moment /iRh = —0.48/1^. This explains the reduced 
average total moment /Z5 = O.S/xb and the very small average Rh moment /XRh = 0.15//S 
found in SP isomer of FeRh4. The present example illustrates the subtle competition between 
cluster structure and magnetism in 3(i-4(i nanoalloys. 

Further increase in the Fe content does not change the topology of the optimal structure. 
Moreover, we start to see that for nearly equal concentrations of Fe and Rh (i.e., Fe2Rh3 and 
Fe3Rh2) the low-lying isomers are the result of changes on the chemical order, i.e., changes in 
the distribution of the Fe and Rh atoms within the cluster, rather than the result of changes 
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TABLE IV. (Color online) Structural, electronic and magnetic properties of FeRh pentamers as in 
Table IE 



Cluster 


Struct. 


Eb 


AEm 










1^0 


Fes 




2.51 


1.00 


2.41 


3.20 


2.93 




227 


Fe4Rli 




2.76 


1.01 


2.30 
2.47 


3.00 


3.09 


1.06 


243 


FesRhs 




2.96 


0.92 


2.37 
2.39 


2.40 


3.13 


0.75 


244 


FeaRhg 




3.06 


0.55 


2.35 
2.71 


2.20 


3.36 


1.08 


260 


FeiRh4 




3.01 


0.33 


2.39 
2.51 


1.20 


3.31 


0.57 


251 


Rhs 




3.03 


0.70 


2.48 


1.00 




0.95 


113 



in the cluster topology. The most stable configuration corresponds to the case where the 3 
Rh atoms (in Fe2Rh3) or the 3 Fe atoms (in Fe3Rh2) are all NNs of each other (see Table 
IV). This is understandable from a single-particle perspective, since the band energy is lower 
when orbitals having nearly the same energy levels are hybridized. In addition, the most 
stable configurations maximize first the number of FeRh NN pairs, followed by the number of 
RhRh pairs .'^Finally, in the Fe-rich limit, for example in Fe4Rh, the lowest-energy structure 
remains a TBP but the closest isomer corresponds to the SP, which has a different topology, 
rather than a different position of the Rh atom in the TBP. 

The trends in the magnetic properties are dominated by the Fe content. As for smaller 
clusters the average magnetic moment per atom /Z^ increases monotonously with increasing 
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Fe concentration. This holds for all optimal structures and in most of the first excited 
isomers. In fact the latter show in general the same /Z^ as the optimal structure. The only 
exception is FeRh4, which is also the only case where an antiparallel alignment of Rh local 
moments is found. In all other investigated cases the magnetic order was found to be FM- 
like. The local Fe moments show the usual enhancement with respect to pure Fcat, due to 
an increase in the number of Fe (i-holes. This effect is stronger for Rh-rich clusters, since the 



larger the number of Rh atoms is, the stronger is the FeRh charge transfer (see Table IV). 
In contrast, the substitution of Rh by Fe does not always enhances the Rh local moments, 
as we observed systematically for smaller sizes. Finally, it is interesting to observe that the 
different chemical orders found in the low lying isomers of Fe2Rh3 and Fe3Rh2 correspond 
to different local magnetic moments. The environment dependence of /Xq, follows in general 
the well-known trend of higher spin polarization at the lowest coordinated sites. 



E. FeRh hexamers 

In Table |V] the results for FeRh hexamers are summarized. For each composition all 
possible cluster topologies (63 different graphd^SED) ^nci all non-equivalent distributions of 
Fe and Rh atoms were taken into account as initial guess for the ab initio optimization of the 
cluster geometry. As in previous cases, all relevant values of the total magnetic moment 23^ 
are scanned. Despite the diversity of starting topologies most low coordinated structures 
relax into compact ones in the course of the unconstrained relaxations. In the end, the 
square bipyramid (SBP), in general somehow slightly distorted, yields the lowest energy 
regardless of composition. 

The binding energy per atom Eb shows a similar composition dependence as for pen- 
tamers. For Fe-rich clusters Eb increases steadily with increasing Rh content, by about 
0.2 eV each time a Rh replaces an Fe (see Table [V|) . Qualitatively, this confirms that the 
bonding between Fe and Rh is stronger than between Fe atoms. However, for nearly equal 
concentrations and in the Rh-rich clusters (FcmRhe-m with m < 2) Eb becomes almost 
independent of m. This seems to be the result of a compensation of bonding and magnetic 
contributions. In fact, on the one side, the magnetic energy AEm continues to decrease 
with increasing Rh content, by about 0.1-0.2 eV per Rh substitution, even for high Rh 
content. And on the other side, this is compensated by an increase of the bonding energy 
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TABLE V. (Color online) Structural, electronic and magnetic properties of FeRh hexamers as in 
Table IE 
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Struct. 
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1.03 


2.38 


3.33 


2.95 
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Fe5Rh 




2.94 


1.07 


2.38 
2.46 


3.16 


3.10 


1.20 


210 


Fe4Rli2 




3.14 


0.85 


2.38 
2.45 


3.00 


3.28 


1.22 


214 


FesRhg 




3.21 


0.56 


2.62 
2.38 
2.61 


2.50 


3.32 


1.14 


239 


Fe2Rli4 




3.26 


0.38 


2.46 
2.51 


2.33 


3.38 


1.39 


207 


FeRhs 




3.24 


0.28 


2.44 
2.54 


1.83 


3.37 


1.33 


192 


Rhe 




3.20 


0.19 


2.54 


1.00 




0.91 


188 



with increasing number of Rh atoms. 

In the case of Rhg an octahedron with an average moment /ig = IfiB and average bond 
length d = 2.54 A yields the lowest energy. The first isomer, a trigonal biprism (TBP), 
lies only 28 meV above the optimum, showing a somewhat shorter average bond length 
d = 2.46 A and a higher average moment JIq = 1.67 fx b- These results are consistent with 
previous DFT calculations.EH A single Fe substitution enhances the average moment to 
Jlq = 1.83/iB but does not change the topology of the optimal FeRhs. The RhRh distances 
remain essentially unchanged and the FeRh distances are somewhat shorter. The important 
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increase in the ground-state spin polarization (S/i^ in all) is not only due to the larger Fe 
moment (/ipe = 3.37/iB in the PAW sphere) but also results from the enhancement of the 
local Rh moments (/iRh = 1.33/1^, see Table [V|). The first isomer of FeRhs corresponds to 
a distorted trigonal prism with significantly contracted FeRh and RhRh bond lengths. In 
Fe2Rh4 the Fe atoms occupy the opposite apex positions of the octahedron. In this way each 
Fe is four-fold coordinated with all Rh atoms. Heterogeneous bonds are favored over FeFe 
ones. The local Fe moments in Fe2Rh4 are the largest among all hexamers: /ipe = 3.38/xb, 
slightly beyond the value found in FeRhs. This corresponds to a large number of d holes. In 
addition, particularly important spin polarizations are induced at the neighboring Rh atoms 
(/^Rh = 1.39/iB). The first isomer of Fe2Rh4 corresponds to a capped trigonal bipyramid 
(CTBP) having a short FeFe NN bond. This structure lies only 0.11 eV higher in energy 
and has the same total moment as the optimal geometry. Replacing a further Fe atom yields 
FeaRhs, whose optimal structure is an octahedron. Here we find two isosceles open Fes and 
Rhs triangles that form a 7r/2 angle with respect to each other (see Table [V|) . Out of the 
12 NN pairs in the FesRhs octahedron, 8 are FeRh and only 4 are homogeneous (2 FeFe 
and 2 RhRh). The local Fe magnetic moments are similar to the other clusters but the 
Rh moments are somewhat smaller in average (/iRh = l-14/i_B). The first excited isomer of 
FeaRhs is a CTBP that lies 25 meV per atom above the ground state. The lowest-energy 
structure found for Fe4Rh2 is an octahedron, while a distorted CTBP is an isomer lying 
0.14 eV per atom above. In the former the Rh atoms are far apart occupying the apical 
positions, whereas in the latter they are NNs. The situation is thus similar to what we find 
in Fe3Rh2. For low Rh or Fe concentrations the atoms are distributed in order to favor the 
FeRh bonds rather than homogeneous NN pairs between the atoms in the minority. The 
octahedron and a distorted CTBP remain the two most stable structures as one further 
reduces the Rh content (see Table for FegRh and Fcg). 

Concerning the magnetic properties one observes qualitatively similar trends as in the 
smaller clusters. The average magnetic moment per atom /Ijy increases monotonously with 
Fe content. Accordingly, the energy gain AE^ associated to magnetism also increases with 
the number of Fe atoms. There are in general very little differences in JIj^ between the 
optimal structure and the first low-lying isomer. The largest part of the spin polarization 
(about 90%) can be traced back to the local d magnetic moments with the PAW sphere of 
the atoms. As expected, the s and p spin polarizations are almost negligible in comparison 
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to the (i-orbital contributions. A significant increase of tlie Fe moments is observed upon Rh 
doping, wliicfi result from tlie larger number of available Fe d holes and the low coordination 
number. Moreover, the Rh moments in FemRh„ are stabilized by the proximity of the Fe 
atoms. In the alloy hexamers the values of /iRh are larger than in the pure Rhg. However, this 
is not a general trend, since the magnetic moments in small Rh„ are often quite important 
due to the extremely reduced coordination numbers. 



F. Exploring heptamers and octamers 

For FcmRh^ clusters having m + n = N > 7 we did not attempt to perform a systematic 
sampling of initial topologies for further unconstrained structural relaxation, as was done 
for the smaller sizes. Instead of aiming at a true global optimization, only a few compact 
and open starting structures are considered. For N = 7, the topologies include bicapped 
trigonal bipyramid (BCTBP), capped octahedra (CO), and pentagonal bipyramid (PBP), 
while for N = 8, they are the tricapped trigonal bipyramid (TCTBP), bicapped octahedra 
(BCO), capped pentagonal bipyramid (CPBP), and cube (C). This choice is motivated 
by previous results for pure clusters and by the trend to compact geometries observed for 
smaller sizes < 6. Although far from exhaustive, the considered geometries allow to 
explore various relevant growth patterns with a reasonable computational effort. Certainly, 
a more complete study would be necessary in order to draw definitive conclusions about the 
optimal topologies. For each composition, all possible distributions of the Fe and Rh atoms 
within the cluster, as well as all relevant values of the total magnetization 5*^ are taken into 
account (from the non-magnetic state to saturation). Therefore, the trends on the interplay 
between chemical order and magnetic behavior remain rigorous within the framework of the 
sampled topologies. 



The results for N = 7 are summarized in Table VI As in smaller FemRhn the binding 
energy per atom increases first with increasing Rh content and becomes essentially indepen- 
dent of composition in the Rh-rich limit (5 < ri < 7). For pure Rh7 the PBP is the most 
stable structure among the considered starting geometries. This result is consistent with 
some earlier DFT studies.'^ However it contrasts with the calculations by Wang et a/.,'^who 
used the GGA functional of Ref. |3l| (PW91) and obtained a capped octahedra, and with the 
calculations of Bae et al.,^"^ who found a prism plus an atom on a square face. According to 

19 



TABLE VI. (Color online) Structural, electronic and magnetic properties of FeRh heptamers as 
obtained from a restricted representative sampling of cluster topologies (see text). 
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Fe4Rh3 




3.36 
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Fe3Rh4 
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0.57 
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1.29 


220 


Rh7 




3.33 


0.22 


2.61 


1.86 




1.62 


203 



our results, these structures are, respectively 20 and 4 meV per atom higher in energy than 
the PBP. In the case of a prism plus an atom on the square face, the energy difference with 
the ground state seems too small to be able to draw definitive conclusions. 

FeRh heptamers with high Rh concentrations also favor a PBP topology. In FeRhg the 
Fe atom occupies an apex site, while in Fe2Rh5 and Fe3Rh4 the Fe atoms belong to the 
pentagonal ring. Notice that the distances between the two apex atoms in Rhy and between 
the Fe and Rh apex atoms in FeRhg are relatively short. In Fe2Rh5 and Fe3Rh4 the Fe 
atoms are as far as possible from each other and the distance between the Rh apex atoms 
is larger. This is consistent with the previously discussed trend to favor the stronger FeRh 
bonds. For example, the energy involved in changing the position of the Fe atom in FeRhg 
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from the apex (7 FeRh bonds) to the pentagonal ring (5 FeRh bonds) is 0.016 eV per atom. 

As the Fe content increases, the topology of FemRh„ changes. In fact Fe4Rh3 corresponds 
to a CO, while for m > 5 the configuration yielding the lowest energy can be regarded as a 



strongly distorted PBP (see Table VI). Already in Fe4Rh3, but also in Fe5Rh2, one observes 
a tendency of the Fe atoms to group in subclusters, bringing the Rh atoms to outer positions, 
so that the number of FeRh bonds is largest. Concerning the shape of the Fe rich heptamers, 
one observes important deformations of the pentagonal bipyramid (D5/1 symmetry) which 
are similar to the distortions found in pure Fey.'™' While the precise origin of the symmetry 
lowering is difficult to establish in the alloys, it is reasonable to expect that it is similar to 
the case of pure Fey. According to Ref. [HI the deformations found in Fe? are due to the 
presence of degenerate electronic states in the undistorted PBP structure. In order to verify 
this hypothesis we have analyzed the Kohn-Sham spectrum of the symmetric structure {D^h 
symmetry) and found that it is highly degenerate at ep- In contrast the spectrum of the 
distorted structure has a band gap about 0.4 eV at ep- This suggests that the distortions 
in Fe rich FeRh clusters can be interpreted as a Jahn- Teller effect. 



In Table VII results for FeRh octamers are reported. The general trends concerning the 
composition dependence of the binding energy, chemical order, as well as the average and 
local magnetic moments are very similar to smaller clusters. The most stable structure 
that we obtain for Feg is a BCO having Eb = 3.03 eV, an average magnetic moment 



/ig = 3jjB, and a relatively short average bond- length d = 2.42 A (see Table VII). A similar 
structure is also found in previous spin-polarized LDA calculations, where Eb = 4.12 eV 
and Jig = 3jjB were obtained.!^ We have repeated these calculation for the BCO structure 
with our computational parameters and atomic reference energies and found Eb = 3.51 eV. 
The discrepancies between LDA and GGA results reflect the importance of exchange and 
correlation to the binding energy. In the other extreme, for pure Rhg, the structure that 
we find with the considered starting topologies is a regular cube having Eb = 3.59 eV, an 
average magnetic moment /Ig = 1.5/2^, and all bond lengths equal to 2.40 A. These results 
are in good agreement with previous calculations by Bae et al..^^ It is interesting to observe 
that the substitution of a single Rh atom by Fe in FeRhy results in a compact topology, 
which is more stable than the relatively open (relaxed) cube-like structures derived from 
pure Rhg. The same trend holds for higher Fe content (i.e., FemRhg_m with m > 1). The 
dominant structure for non-vanishing Fe content is a BCO with slight distortions. Only for 
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TABLE VII. (Color online) Structural, electronic and magnetic properties of FeRh octamers as 
obtained from a restricted representative sampling of cluster topologies (see text). 
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FesRhs we find a different topology, namely, a distorted CPBP. The typical isomerization 
energies between the BCO and the TCTBP are AE'iso = 10-30 meV per atom. The average 
magnetic moments in the lowest lying isomers are either the same or very similar. 

The magnetic properties of heptamers and octamers follow qualitatively the behavior 
observed in smaller clusters. In most cases the average magnetic moment per atom JIj^ 
and the magnetic energy AEm increase with Fe concentration. The only exception is the 
pure Rh heptamer, for which /I^ is somewhat larger than in FeRhg. This is not due to 
AF-like coupling between Fe impurity moment and the remaining Rh atoms but rather to a 
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reduction of the Rh local moments in FeRhg (/uf^ — 1.61-1.63/iB in Rhy, while fif^ ~ 1.25- 
in FeRhg). Remarkably, the Rh local moments in Rhy are the largest among all the 
heptamers. They amount to 87% of the total moment, which stresses the importance of the 
local (i-electron contributions. Also in Rhg one finds quite large local moments, which are 
actually larger than the Rh moments in most Fe doped clusters. This shows that for these 
sizes the Fe atoms do not necessarily increase the Rh moments by simple proximity effects 



(see Tables VI and VII). Nevertheless, a different behavior is expected for larger A^, where 
pure Rh clusters are no longer magnetic on their own. The local Fe moments are strongly 
enhanced with respect to pure FeAr (/xpe — 2.8/iB in Fey or Feg) reaching values up to 3.36/2^, 
particularly when the Fe atoms are in a Rh rich environment. As in the smaller clusters, 
this is a consequence of a charge transfer from the Fe to the Rh atoms, which increases 
the number of polarizable Fe d-holes. Notice that some kind of interaction between the Fe 
atoms seems to favor this effect, since the largest /ipe are found for clusters having 2 or 3 
Fe atoms rather than for the single Fe impurity. Large Fe moments are also found in bulk 
FeRh alloys.li^ 

To conclude this section it is interesting to compare the cluster results with available 
experiments and calculations for macroscopic alloys Band structure calculations for 
the periodic Feo.sRho.s alloy having a CsCl structure yield an antiferromagnetic (AF) ground 
state, which is more stable than the ferromagnetic solution.'^ This is qualitatively in agree- 
ment with experiments showing AF order when the Rh concentration is above or equal to 
50%^ In contrast our results for small clusters show a FM-like order for all Rh concentra- 
tions, even for the pure Rh clusters. This is a consequence of the reduction of local coor- 
dination number and the associated effective rf-band narrowing, which renders the Stoner 
criterion far easier to satisfy, and which tends to stabilize the high-spin states with respect to 
the low-spin AF states. In fact, even in the bulk calculations on FeRh, the energies of the AF 
and FM states are not very different, and a coexistence of both solutions is found over a wide 
range of volumes.l^ Moreover, experiment shows an AF to FM transition with increasing 
temperature, which is accompanied by an enhanced thermal expansion.!^ Recent ab initio 
calculations have revealed the importance of competing FM and AF exchange interactions 
in stoichiometric a-FeRh.'^ Moreover, neutron diffraction experiment^ on Fei_a;Rha; for 
0.35 < X < 0.5 and calculationJ^^ for x = 0.5 show that the Fe moments /ipe are significantly 
enhanced with respect to /xpe in pure a-Fe, particularly in the FM state where it reaches 
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values of about 3.2/xb.'2E3 These bulk results are remarkably similar to the trends found 
in FemRh„ clusters over a wide range of compositions. As in the clusters, the induced Rh 
moments /XRh play an important role in the stability of the FM phase. Bulk experiment^ 
on Fei-xBhx yield /XRh — for 0.35 < x < 0.5 which is comparable to, though somewhat 
smaller than the present cluster results. 

IV. SIZE AND COMPOSITION DEPENDENCE 

The main purpose of this section is to focus on the dependence of the electronic and mag- 
netic properties of FemR-h„ clusters as a function of size and composition. In the following 
we present and discuss results for the binding energy, average and local spin moments, and 
electronic densities of states for N = m + n < 8. 

A. Binding energy and magnetic moments 

In Fig. [T] the binding energy per atom Eb is given as a function of the number of Fe 
atoms m. Besides the expected monotonic increase of Eb with increasing N, an interesting 
concentration dependence is observed. For very small sizes (A^ < 4) Eb is maximal for 
m = 1 or 2, despite the fact that Eb is always larger for pure Rh than pure Fe clusters. 
This indicates that in these cases the bonding resulting from FeRh pairs is stronger than 
RhRh bonds. Only for m > N — 1, when the number of weaker FeFe bonds dominates, 
one observes that Eb decreases with increasing m. For larger sizes (A^ > 5) the strength of 
RhRh and FeRh bonds becomes very similar, so that the maximum in Eb is replaced by a 
range of Fe concentrations x = m/N < 0.5 where Eb depends very weakly on m. 

In Fig. |2] the average magnetic moments Jtj^ of FemRh„ are shown as a function of m 
for A^ < 8. As already discussed in previous sections, JIj^ increases monotonously, with the 
number of Fe atoms. This is an expected consequence of the larger Fe local moments and 
the underlying FM-like magnetic order. The average slope of the curves tends to increase 
with decreasing A^, since the change in concentration per Fe substitution is more important 
the smaller the size is. The typical increase in /I^ per Fe substitution is about {1/N)^b per 
Fe substitution. Notice, moreover, the enhancement of the magnetic moments of the pure 
clusters in particular for FeAr (m = A^), which go well beyond S/^s, the value corresponding 
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FIG. 1. (Color online) Binding energy per atom Eb of FemRhn clusters as a function of the number 
of Fe atoms. The lines connecting the points for each N = m + n are a guide to the eye. 
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FIG. 2. (Color online) Total magnetic moment per atom Jlj^j of FemR.h„ clusters as a function of 
number of Fe atoms. The symbols corresponding to each size are the same as in Fig. [TJ The lines 
connecting the points for each N = m + n are a guide to the eye. 

to a saturated d-hand in the d'^s^ configuration. In contrast, the moments of pure Rlijv are 
far from saturated except for = 2 and 7 (see Fig. |2] for m = 0). In this context it is 
important to recall that a thorough global optimization, for example, by considering a large 
number of initial topologies, could affect the quantitative values of the magnetic moments 
for = 7 and 8. 

The local magnetic moments in the PAW sphere of the Fe and Rh atoms provide further 
insight on the interplay between 3d and Ad magnetism in Fe^Rh^. In Fig. |3] /xpe and finh 
are shown as a function of m for A^ = 6-8. The Fe moments are essentially given by the 
saturated ci-orbital contribution. For pure Fe clusters the actual values of /xpe within the 
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FIG. 3. (Color online) Local magnetic moment [la at the Fe and Rh atoms as a function of the 
number Fe atoms m. 



PAW sphere are somewhat lower than S/i^ due to a partial spill-off of the spin-polarized 
density. Notice that the Fe moments increase as we replace Fe by Rh atoms showing some 
weak oscillations as a function of m. The increase is rather weak for a single Rh impurity 
in FcAT-iRh but becomes stronger reaching a more or less constant value as soon as the 
cluster contains 2 or more Rh (m < — 2, see Fig. [3]). This effect can be traced back to a c? 
electron charge transfer from Fe to Rh which, together with the extremely low coordination 
number, yields a full polarization of the larger number of Fe d holes. On the other side 
the Rh moments are not saturated and therefore are more sensitive to size, structure and 
composition. The values of /iRh are in the range of 1-1. 5/^^ showing some oscillations as a 
function of m. No systematic enhancement of /iRh with increasing Fe content is observed. 
This behavior could be related to charge transfers effects leading to changes in the number 
of Rh d electrons function of m. 

Finally, it is interesting to analyze the role played by magnetism in defining the cluster 
structure by comparing magnetic and non-magnetic calculations. For the smallest FeRh 
clusters (iV = 3 and 4) the magnetic energy l^Em = E{Sz = 0) — E{Sz) gained upon 
magnetization is higher in the first excited isomer than in the most stable structure. This 
implies that the contribution of magnetism to the structural stability is not crucial, since 
the non-magnetic calculations yield the same ordering, at least concerning the two best 
structures. This suggests that for the smallest sizes the kinetic or bonding energy dominates 
the structural stability, which also explains that the two most stable isomers have different 
topologies. The situation changes for large clusters. For > 5 one finds a number of FeRh 
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clusters for which the optimal structure is actually stabilized by magnetism. For example, 
in Fe4Rh, Fe3Rh2, and FeRh4 the energy ordering of the two most stable isomers would be 
reversed if magnetism were neglected. It should be noted that in these cases the structures 
differ only in the chemical order, not in the topology which is a TBP. In the FeRh hexamers 
the energy differences between the low-lying isomers are more important and only in one 
case, Fe4Rh2, magnetism appears to be crucial for stabilizing the actual optimal structure. 
A similar strong interplay between structure, chemical order and magnetism is expected for 
larger FeRh clusters. 



B. Electronic structure 

In the previous sections the structure and spin moments of FeRh alloy clusters have been 
discussed as a function of 3d/ Ad concentration. Although these properties are intimately 
related to the size and composition dependence of electronic structure, it is in general very 
difficult to achieve a physical transparent correlation between global and microscopic behav- 
iors. Nevertheless, it is very interesting to analyze, at least for some representative examples, 
how the electronic structure depends on the composition of magnetic nanoalloys. To this 
aim we report in Fig. |4]the spin-polarized li-electron density of states (DOS) of representa- 



tive FeRh octamers having the relaxed structures illustrated in Table VII Results for pure 
Feg and Rhg are also shown for the sake of comparison. In all the clusters, the dominant 
peaks in the relevant energy range near ep correspond either to the Fe-3d or to the Rh-4(i 
states. The valence spectrum is largely dominated by these ci-electron contributions. In fact 
the total DOS and the (i-projected DOS are difficult to tell apart. 

First of all, let us consider the DOS of the pure clusters. Our results for Rhg with a 
cube structure are similar to those of previous studies.!^ They show the dominant (i-electron 
contribution near ep, with the characteristic ferromagnetic exchange splitting between the 
minority and majority spin states. In Fig. |4] we also included the DOS for Rhg with a 
BCO structure, since it allows us to illustrate the differences in the electronic structure 
of compact and open geometries. Moreover, the DOS of pure Rhg with BCO structure 
is very useful in order to demonstrate the dependence of DOS on Fe content, since the 
structures of Fem,Rhg_m with m > 1 are similar to the BCO. Both Feg and Rhg show 
relatively narrow d-bands which dominate the single-particle energy spectrum in the range 
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FIG. 4. (Color online) Electronic density of states (DOS) of FeRh octamers. Results are given 
for the total (solid), the Fe-projected (dotted), and the Rh-projected (dashed) d-electron DOS. 
Positive (negative) values correspond to majority (minority) spin. A Lorentzian width A = 0.02 eV 
has been used to broaden the discrete energy levels. The corresponding structures are illustrated 
in Table IVlIl 

— 5eV < s — sp < 3eV. The spin polarization of the DOS clearly reflects the ferromagnetic 
order in the cluster. Putting aside the exchange splitting, the peak structure in the up and 
down DOS Ptj{e) are comparable. There are even qualitative similarities between the two 
elements. However, looking in more detail, one observes that the effective (i-band width in 
Feg (about 4 eV) is somewhat smaller than in Rhg (about 5 eV). Moreover, in Rhg the DOS 
at eiT' is non- vanishing for both spin directions and the finite-size gaps are very small (see 
Fig. |4]). In contrast, the majority ci-DOS is fully occupied in Feg, with the highest majority 
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state lying about 0.5 eV below sp- In addition there is an appreciable gap (about 0.1 eV) in 
the corresponding minority spectrum. These qualitative differences are of course consistent 
with the fact that Feg is a strong ferromagnet with saturated moments, while Rhg should 
be regarded weak unsaturated ferromagnet. 

The trends as a function of concentration reflect the crossover between the previous 
contrasting behaviors. For low Fe concentration (e.g., Fe2Rh6) we still find states with both 
spin directions close to Sp- The magnetic moments are not saturated, although the Fermi 
energy tends to approach the top of the majority band. Moreover, the majority-spin states 
close to Bp have dominantly Rh character. Small Fe doping does not reduce the d-band 
width significantly. Notice the rather important change in the shape of the DOS in Fe2Rh6 
as compared to the DOS in Rhg. This is a consequence of the change in topology from cubic 
to bicapped octahedron (BCO). 

For equal concentrations (Fe4Rh4) the first signs of rf-band narrowing and enhanced 
exchange splitting start to become apparent. The spin- up states (majority band) which in 
Fe2Rhe contribute to the DOS a,t sp now move to lower energies (0.3 eV below Sp) so that 
the majority band is saturated. Only spin-down (minority) states are found around Sp, 
although there is a significant gap in pi{e) (see Fig. |4]). In the majority band Rh dominates 
over Fe at the higher energies (closer to Sp), while Fe dominates in the bottom of the band. 
In the minority band the participation of Rh (Fe) is stronger (weaker) below Sp and weaker 
(stronger) above ep. This is consistent with the fact that the Rh local moments are smaller 
than the Fe moments. 

Finally, in the Fe rich limit (e.g., Fe6Rh2), the majority-band width becomes as narrow as 
in Feg, while the minority band is still comparable to Rhg. The exchange splitting is large, 
the majority band saturated and only minority states are found close to Sp. As in Feg, Pi{e) 
shows a clear gap aX Sp (see Fig. |4]). However, the Rh contribution to the minority states 
below Bp remains above average despite the relative small Rh content. The Fe contribution 
largely dominates the unoccupied minority-spin DOS, in agreement with the larger local Fe 
moments. 
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V. DISCUSSION 



The structural, electronic and magnetic properties of small Fe^Rh^ clusters having N — 
m + n < 8 atoms have been investigated systematically in the framework of a generalized 
gradient approximation to density-functional theory. For very small sizes (A^ < 4 atoms) the 
binding energy Eb shows a non-monotonous dependence on concentration, which implies 
that the FeRh bonds are stronger than the homogeneous ones. However, for larger sizes 
the FeRh and RhRh bond strengths become comparable, so that Eb depends weakly on 
concentration for high Rh content. 

The magnetic order of the clusters having the most stable structures is found to be FM- 
like. Moreover, the average magnetic moment per atom /Z^ increases monotonously, which 
is almost linear over a wide range of concentration with Fe content. Consequently, the 
energy gain AEj^ associated to magnetism also increases with the number of Fe atoms. The 
largest part of the spin polarization (about 90%) can be traced back to the local d magnetic 
moments within the PAW sphere of the atoms. The s and p spin polarizations are almost 
negligible in general. A remarkable enhancement of the local Fe moments is observed as 
result of Rh doping. This is a consequence of the increase in the number of Fe d holes, due 
to charge transfer from Fe to Rh, combined with the extremely reduced local coordination. 
The Rh local moments are important already in the pure clusters (A^ < 8). Therefore, they 
are not significantly enhanced by Fe doping. However, the overall stability of magnetism, as 
measured by the total energy gained by the onset of spin polarization, is found to increase 
with increasing Fe content. 

A further interesting aspect, particularly for future studies is the role of spin-orbit (SO) 
interactions on the magnetism of nanoalloys. We have performed some representative calcu- 
lations by taking into account spin-orbit coupling (SOC) in order to explore their effect on 
the ground-state structure, chemical order and spin moments. For example in Feg, FcaRha 
and Rhe we find that the changes in the ground-state energy resulting from SO interaction 
are typically of the order of 0.2 cV for the whole cluster. This is often comparable to or larger 
than the energy differences between the low-lying isomers. However, the SO energies are 
very similar for different structures, so that the ground-state structures remain essentially 
the same as in the scalar relativistic (SR) calculations. The changes in the bond lengths and 
in the spin moments resulting from SOC are also very small (e.g., j/Z^*^*-' — ~ 0.01/xb 
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and \dfp^ — dfj^l ~ O.OOlA in FesRhs). As a result, the conclusions drawn from our SR 
calculations on the relative stability and local spin moments seem to be unaffected by the 
spin-orbit contributions. 

FeRh clusters are expected to develop a variety of further interesting behaviors, which still 
remain to be explored. For instance, larger cluster should show a more complex dependence 
of the magnetic order as a function of concentration. In particular for large Rh content one 
should observe a transition from FM-like to AF-like order with increasing cluster size, in 
agreement with the AF phase found in solids for more than 50% Rh concentration. Moreover, 
the metamagnetic transition observed in bulk FeRh alloys also puts forward the possibility 
of similar interesting phenomena in nanoalloys as a function of temperature. Finally, the 
contributions of orbital magnetism and magnetic anisotropy deserve to be explored in detail 
as a function of composition and chemical order, even for the smallest sizes, particularly 
because of their implications for potential applications.'^ 
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